Spontaneous decay dynamics in atomically doped carbon nanotubes 
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We report a strictly non-exponential spontaneous decay dynamics of an excited two-level atom 
placed inside or at diflterent distances outside a carbon nanotube (CN). This is the result of strong 
non-Markovian memory effects arising from the rapid variation of the photonic density of states with 
frequency near the CN. The system exhibits vacuum-field Rabi oscillations, a principal signature 
of strong atom-vacuum-field coupling, when the atom is close enough to the nanotube surface and 
the atomic transition frequency is in the vicinity of the resonance of the photonic density of states. 
Caused by decreasing the atom-field coupling strength, the non-exponential decay dynamics gives 
place to the exponential one if the atom moves away from the CN surface. Thus, atom-field coupling 
and the character of the spontaneous decay dynamics, respectively, may be controlled by changing 
the distance between the atom and CN surface by means of a proper preparation of atomically 
doped CNs. This opens routes for new challenging nanophotonics applications of atomically doped 
CN systems as various sources of coherent light emitted by dopant atoms. 

PACS numbers; 61.46.-|-w, 73.22.-f, 73.63.Fg, 78.67.Ch 



I. INTRODUCTION 

It has long been recognized that the spontaneous emis- 
sion rate of an excited atom is not an immutable prop- 
erty, but that it can be modified by the atomic envi- 
ronment. Generally called the Purcell effect the phe- 
nomenon is qualitatively explained by the fact that the 
local environment modifies the strength and distribution 
of the vacuum electromagnetic modes with which the 
atom can interact, resulting indirectly in the alteration 
of atomic spontaneous emission properties. 

The Purcell effect took on special significance recently 
in view of rapid progress in physics of nanostructures. 
Here, the control of spontaneous emission has been pre- 
dicted to have a lot of useful applications, ranging from 
the improvement of existing devices (lasers, light emit- 
ting diodes) to such nontrivial functions as the emis- 
sion of nonclassical states of light 2]. In particular, the 
enhancement of the spontaneous emission rate can be 
the first step towards the realization of a thresholdless 
laser or a single photon source 0]. The possibility to 
control atomic spontaneous emission was shown theoret- 
ically for microcavities and microspheres H, |^ , optical 
fibers photonic crystals 0], semiconductor quantum 
dots lUJ . Recent technological progress in the fabrication 
of low-dimensional nanostructures has enabled the exper- 
imental investigation of spontaneous emission for micro- 
cavities photonic crystals (Q, quantum dots jlSj . 

Carbon nanotubes (CNs) are graphene sheets roUed- 
up into cylinders of approximately one nanometer in di- 
ameters. Extensive work carried out worldwide in re- 
cent years has revealed the intriguing physical proper- 
ties of these novel molecular scale wires hj hj. Nan- 
otubes have been shown to be useful for miniaturized 
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electronic, mechanical, electromechanical, chemical and 
scanning probe devices and materials for macroscopic 
composites [16j . Important is that their intrinsic prop- 
erties may be substantially modified in a controllable 
way by doping with extrinsic impurity atoms, molecules 
and compounds |17| . This opens routes for new chal- 
lenging nanophotonics applications of atomically doped 
CN systems as various sources of coherent light emitted 
by dopant atoms. Recent successful experiments on en- 
capsulation of single atoms into single-wall carbon nan- 
otubes lia and their intercalation into single-wall CN 
bundles |l7l 0| stimulate an analysis of atomic sponta- 
neous emission in such systems as a first step towards 
their nanophotonics applications. 

Typically, there may be two qualitatively different 
regimes of interaction of an atomic excited state with 
a vacuum electromagnetic field in the vicinity of the CN. 
They are the weak coupling regime and the strong cou- 
pling regime ^2d\. The former is characterized by the 
monotonous exponential decay dynamics of the upper 
atomic state with the decay rate altered compared with 
the free-space value. The latter is, in contrast, strictly 
non-exponential and is characterized by reversible Rabi 
oscillations where the energy of the initially excited atom 
is periodically exchanged between the atom and the field. 
In the present paper, we develop the quantum theory of 
the spontaneous decay of an excited two-level atom near 
a CN and derive the evolution equation of the upper state 
of the system. By solving it numerically, we demonstrate 
the strictly non-exponential spontaneous decay dynamics 
in the case where the atom is close enough to the CN sur- 
face. In certain cases, the system exhibits vacuum-field 
Rabi oscillations — a result already detected for quasi- 
2D excitonic and intersubband electronic transitions in 
semiconductor quantum microcavities pH l22j | and never 
reported for atomically doped CNs so far. 

The rest of the paper is arranged as follows. Section lU 
presents a theoretical model we use to derive the evolu- 
tion equation of the upper state population probability 
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FIG. 1: (Color online) The geometry of the problem. 



amplitude of the composed quantum system " a two-level 
atom interacting with a quantized CN-modified vacuum 
radiation field" . In describing atom-field interaction we 
follow the original line of Refs. adopting their elec- 

tromagnetic field quantization scheme for a particular 
case of the field near an infinitely long achiral single- 
wall CN. The carbon nanotube is considered as an in- 
finitely thin anisotropically conducting cylinder. Its sur- 
face conductivity is represented in terms of the 7r-electron 
dispersion law obtained in the tight-binding approxima- 
tion with allowance made for azimuthal electron momen- 
tum quan tization and axial electron momentum relax- 
ation [231 . Only the axial conductivity is taken into ac- 
count and the azimuthal one, bein g st ron gly supp ressed 
by transverse depolarization fields l25ll26l l27j . is ne- 
glected. In Section ITTll the time dynamics of the upper 
state population probability is analyzed qualitatively in 
terms of two different approximations admitting analyt- 
ical solutions of the evolution equation derived in Sec- 
tion ^] They are the Markovian approximation and the 
single-resonance approximation of the density of photonic 
states in the vicinity of the CN. Section Hvl presents and 
discusses the results of the numerical solution of the evo- 
lution equation for various particular cases where the 
atom is placed in the center and near the wall inside, 
and at different distances outside achiral CNs of differ- 
ent radii. A summary and conclusions of the work are 
given in Section Ivl 



II. THE MODEL 

The quantum theory of the spontaneous decay of ex- 
cited atomic states in the presence of the CN involves 
an electromagnetic field quantization procedure. Such a 
procedure faces difficulties similar to those in quantum 
optics of 3D Kramers-Kronig dielectric media where the 
canonical quantization scheme commonly used does not 
work since, because of absorption, corresponding opera- 
tor Maxwell equations become non-Hermitian [28j. As a 



result, their solutions cannot be expanded in power or- 
thogonal modes and the concept of modes itself becomes 
more subtle. We, therefore, use an alternative quantiza- 
tion scheme developed in Refs. W^^, where the Fourier- 
images of electric and magnetic fields are considered as 
quantum mechanical observables of corresponding elec- 
tric and magnetic field operators. The latter ones satisfy 
the Fourier-domain operator Maxwell equations modified 
by the presence of a so-called operator noise current den- 
sity J(r,ti;) written in terms of a 3D vector bosonic field 
operator f(r, w) and a medium dielectric tensor e(r, w) 
(supposed to be diagonal) as 

= ^^/nim€u{r,uj) f^{r,uj), i = 1,2,3. (I) 

ZTT 

This operator is responsible for correct commutation re- 
lations of the electric and magnetic field operators in the 
presence of medium-induced absorbtion. In this formal- 
ism, the electric and magnetic field operators are ex- 
pressed in terms of a continuum set of the 3D vector 
bosonic fields f(r,a;) by means of the convolution over r 
of the current with the classical electromagnetic field 
Green tensor of the system. The bosonic field operators 
f^(r,a;) and f(r, w) create and annihilate single-quantum 
electromagnetic medium excitations. They are defined 
by their commutation relations and play the role of the 
fundamental dynamical variables in terms of which the 
Hamiltonian of the composed system "electromagnetic 
field + dissipative medium" is written in a standard sec- 
ondly quantized form as 

H ^ j dv j^dujhuji\v,uj)-i{v,Lo). (2) 

Consider a two-level atom positioned at an arbitrary 
point ta near the CN. Since the problem has a cylin- 
dric symmetry, it is convenient to assign the orthonor- 
mal cylindric basis {6^,6(^,62} in such a way that = 
TA^r = {?'A,0,0} and is directed along the nanotube 
axis (see Figure^. Let the atom interact with a quan- 
tized electromagnetic field via an electric dipole transi- 
tion of frequency ua- The atomic dipole moment may be 
assumed to be directed along the CN axis, d = dz^z- The 
contribution of the transverse dipole moment orienta- 
tion is suppressed because of strong depolarization of the 
transverse field in an isolated CN [2l|2l|2l|23|. Strong 
transverse depolarization along with transverse electron 
momentum quantization allow one to neglect the az- 
imuthal current and radial polarizability 23], in which 
case the dielectric tensor components e,,. and e^p^ are 
identically equal to unit. The component e^^ is caused 
by the CN longitudinal polarizability and is responsible 
for the axial surface current parallel to the vector. This 
current may be represented in terms of the ID bosonic 
field operators by analogy with Eq. JQ. Indeed, taking 
into account the dimensionality conservation in passing 
from bulk to a monolayer in Eq. |2Jl and using a simple 
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Drude relation [2l| 



-ILO- 



(3) 



where R = {i?cri! (j^i Z} is the radius- vector of an arbitrary 
point of the CN surface, Ren is the radius of the CN, 
CTzz(R, w) is the CN surface axial conductivity per unit 
length, 5* is the area of a single nanotube, pT is the tubule 
density in a bundle, one immediately has from Eq. 



J(R,tj) = y /(R,w)ez (4) 

with /(R, oj) being the ID bosonic field operator defined 
on the CN surface. The total Hamiltonian of the system 
under consideration is then written in terms of /^(R, cj) 
and /(R, cj) operators as 

■H = Jd'R-J^ dujhuj P (11,0;) fCR,^) + ^hLu ACT z 

-[a^ Ei+\rA)d,+h.c.]. (5) 

Here, the three terms represent the electromagnetic field 
modified by the presence of the CN, the two-level atom 
and their interaction (within the framework of electric 
dipole, and rotating wave approximations 0), respec- 
tively. The Pauli operators az = \u){u\ — \l){l\, a=\l){u\, 
(T^ = |u) (/| describe the two- level atomic subsystem where 
|u) and 1^) are the upper and lower atomic states, respec- 
tively. The operators E(^^(r^) represent the electric field 
the atom interacts with. For an arbitrary r — {r,ip,z}, 
they are defined as follows 

E(r) =E(+)(r)+E(-)(r), 

E(+)(r)=/ E{r,uj)duj, &-\r) ^ [&+\r)]^ . (6) 
Jo 

Here, E(r, uj) satisfies the Fourier-domain Maxwell equa- 
tions 

V X E(r,Lj) = ikn{r,uj), 

V X H(r, uj) = -ifcE(r, uj) + — I(r, u), (7) 

c 

where H(r, uj) stands for the magnetic field operator [de- 
fined by analogy with Eq. 10)], k — uj/c, and 



I(r, w) =y dR (5(r-R) J(R, uj) = 23 {Ren, ^p, z, uj)S{r-Rcn) 

is the exterior operator current density [with J(R, uj) de- 
fined by Eq. 10}] associated with the presence of the CN. 
From Eqs. Q and © in view of Eq. Q, it follows that 

E(r,a;) =i^fc /"dRG(r,R,w)-J(R,w) (9) 



[and H = {ik)~^\7 x E accordingly], where G(r, R, uj) is 
the Green tensor of the classical electromagnetic field in 
the vicinity of the CN. Its components satisfy the equa- 
tion 

^ (Vx Vx -fc2)^^G„z(r,R,u;)-(5(r-R), (10) 



together with the radiation conditions at infinity and the 
boundary conditions on the CN surface. The Hamilto- 
nian ISJ) along with Eqs. l|^- l(Tn|) and 01 is the modifi- 
cation of the Jaynes-Cummings model 28] for an atom 
in the vicinity of a solitary CN. The classical electro- 
magnetic field Green tensor of this system is derived in 
Appendix ^ 

When the atom is initially in the upper state and the 
field subsystem is in vacuum, the time-dependent wave 
function of the whole system can be written as 

\m) = Cu{t)e-^'---/'^'\u)m) 



dr dujCi{r,u;,t)e-'^^-^-^'^'\l)\{l{r,u;)}), (11) 
Jo 

where [{0}) is the vacuum state of the field subsystem, 
]{l(r, uj)}) is its excited state where the field is in a single- 
quantum Fock state, C'u and C'l are the population proba- 
bility amplitudes of the upper state and lower state of the 
whole system, respectively. In view of Eqs. ©, ©, 
and the integral relationship 



Im Gai3{r,r',uj) 



An 
c 



k J dRReazz (R, c^)G„z (r, R, iu)G},,{r', R, lu) 



(which is nothing but a particular case of the gen- 
eral 3D Green tensor integral relationship rigorously 
proven in Ref. with Eq. Q taken into account), the 
Schrodinger equation with the Hamiltonian jSj) and wave 
function Ijlll) yields the following evolution law for the 
population probability amplitude of the upper state of 
the system 



Cuir) 



1+ / K{T-T')Cu{T')dT', (12) 
"'0 



Kir^r')^^!f^l dxx'a^)- 



i(x-XA)(T-T') _ 2 



i {x - xa) 



Here, 



huj 

2^0 



and 



270* 



(13) 



(14) 



are the dimensionless frequency and time, respectively, 
with 7o = 2.7 eV being the carbon nearest neighbor hop- 
ping integral 29] appearing in the CN surface axial con- 
ductivity in Eq. Fq is the rate of the exponential 
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atomic spontaneous decay in vacuum obtained from the 
general expression of the form [sol Isif 



Fix) 



Im Gzz{T^A,'rA,x) 



by substituting 



1 270 a: 
Gttc Ti 



(15) 



(16) 



for the vacuum imaginary Green tensor The func- 
tion ^(x) is the relative density (with respect to vacuum) 
of photonic states near the CN given for ta > Ren (see 
Appendix IB|| by 



^/ N r(a;) ^ 3 ^ ^ 

ro(x) TT 



(17) 



s(i?cn, a;)w(2/)''/p [i;(y)u(-Rcn)a;]iir2[i;(y)u(r^)a;] 
C 1 + s(i?c„, a:)w(i/)2/p[w(2;)M(i?c„)a:]i4:p[w(2;)M(i?c„)2;] ' 

where /p and Kp are the modified cylindric Bessel func- 
tions, v{y) = - 1 , u(r) = 2jor/hc, and s(i?cn,a;) = 
2iau{Rcn)xazz{Rcn,x) with CTzz = 2Trhazz/e'^ being the 
dimensionless conductivity and a — e^/hc = 1/137 repre- 
senting the fine-structure constant. The integration con- 
tour C goes along the real axis of the complex plane 
and envelopes the branch points y = ±1 of the function 
v{y) in the integrand from below and from above, respec- 
tively. For TA < Rem Eq. (|17l) is modified by the simple 
replacement va Ren in the Bessel function arguments 
in the numerator of the integrand. Note the divergence of 
f (x) at rA — Ren 1 i- e. when the atom is located right on 
the CN surface. The point is that the CN dielectric tensor 
longitudinal component ezz [which, according to Eq. 
is responsible for the surface axial conductivity Gzz in 
Eq. (|17() ] is obtained as a result of a standard procedure of 
physical averaging a local electromagnetic field over the 
two spatial directions in the graphene plane js^ . Such av- 
eraging does not assume extrinsic atoms on the graphene 
surface. To take them into consideration the averaging 
procedure must be modified. Thus, the applicability do- 
main of our model is restricted by the condition 



\rA- Rcn \ > a, 



(18) 



where a = 1.42 A is a graphene interatomic distance [2^ . 

Eq. (|12|) is a well-known Volterra integral equation of 
the second kind. In our case, it describes the sponta- 
neous decay dynamics of the excited two-level atom in the 
vicinity of the CN. All the CN parameters that are rele- 
vant for the spontaneous decay are contained in the rela- 
tive density of photonic states H17(l appearing in the ker- 
nel H13|l . The density of photonic states is, in turn, deter- 
mined by the imaginary classical Green tensor of the CN- 
modified electromagnetic field via r(a::) given by Eq. (|15|l . 



III. QUALITATIVE ANALYSIS 

In this Section we will qualitatively analyze the time 
dynamics of the upper state population probability C'u (t) 
in terms of two different approximations admitting an- 
alytical solutions of the evolution problem ifT^ . itT^ . 
They are the Markovian approximation and the single- 
resonance approximation of the relative density £^{x) of 
photonic states in the vicinity of the CN. 



1. Markovian approximation 

In the case where the Markovian approximation is ap- 
plicable, or, in other words, when the atom-field coupling 
strength is weak enough for atomic motion memory ef- 
fects to be insignificant, so that they may be neglected, 
the time-dependent factor in the kernel (|13|l may be re- 
placed by its long-time approximation 
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i{x - xa) 



-Trd{x — Xa) + iV- 



1 



XA 



{V denotes a principal value). Then, in view of Eq. H17|l . 
the kernel becomes 

/N Kr(xA) . . , N 

K{t - t') = j--^ + 2A(a;A) 



with 

A(xa) = 
and Eq. yields 
C„(t) 



470 



hTo{xA) 



dx X 



•c-xa' 



exp 



K1^{xa) 
470 



i/S.{xA) 



(19) 



— the exponential decay dynamics of the [shifted by 
A(a;yi)] upper atomic level with the rate T{xa)- This case 
was analyzed in Ref. [33| . 



2. Single-resonance approximation of the relative density 
of photonic states 

Another approximation that admits an analytical so- 
lution of the evolution problem H12|l. (|13|l is a single- 
resonance approximation. Suppose that at x = Xr the 
photonic density of states ^(x) has a sharp peak of half- 
width-at-half-maximum Sxr- For all x in the vicinity 
of Xj., the shape of ^(x) may then be roughly approxi- 
mated by the Lorentzian function of the form 



ax) 



^(Xr)Sx'^ 



{x — Xr)"^ + 5x1 

The kernel H13|) is then easily calculated analytically to 
give 



K{t - t') 



hViXr) SXr 
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exp[—i{xr — iSxr — xa){t — t') — 1] 



T > T 



Substituting this into Eq. (|12|l and making the differenti- 
ation of both sides of the resulting equation over time, fol- 
lowed by the change of the integration order and one more 
time differentiation, one straightforwardly arrives at a 
second order ordinary homogeneous differential equation 

(7„(t) + liXr - iSXr - XA)CuiT) + (^1/2)^ ^(t) = 0, 



where — y^26xrKr{xr)/'2,jQ , with the solution given 
for XA ~ Xr by 



Cuir) 




exp 



- ^Sxr— 



1-- 



dXr 



^5x^ - 



exp 



(^Sxr + \/Sa 



(20) 

This solution is approximately valid for those atomic 
transition frequencies xa which are located in the vicin- 
ity of the photonic density-of-states resonances what- 
ever the atom-field coupling strength is. In particular, 
if {^/Sxr)"^ <C 1, Eq. (pn|l yields the exponential de- 
cay of the upper atomic state population probability 
|C'ti(''')P with the rate hr{xr)/2'yQ in full agreement with 
Eq. 119(1 obtained within the Markovian approximation 
for weak atom-field coupling. In the opposite case, when 
(ri/fer)^^!, one has 



\Cu{r)\' 



e-'^"^ cos^ 



and the decay of the upper atomic state population prob- 
ability proceeds via damped Rabi oscillations. This is the 
principal signature of strong atom-field coupling which is 
beyond the Markovian approximation. Expressing T{xr) 
in 17 in terms of £,{xr) by means of Eq. p7|l and using 
the approximation TQ{x)^a'^2jQx/h valid for hydrogen- 
like atoms 35] , one may conveniently rewrite the strong 
atom-field coupling condition in the form 



2a'^x.. 



6Xr 



» 1, 



(21) 



from which it follows that the strong coupling regime is 
fostered by high and narrow resonances in the relative 
density of photonic states. 



IV. NUMERICAL RESULTS AND DISCUSSION 

To get beyond the Markovian and single-peak approx- 
imations, we have solved Eqs. 112(1 and ((13(1 numeri- 
cally. The exact time evolution of the upper state pop- 
ulation probability |C„(t)P was obtained for the atom 
placed [in a way that Eq. ((18() was always satisfied] in the 
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FIG. 2: (Color online) Relative density of photonic states (a) 
and upper-level spontaneous decay dynamics (b) for the atom 
in the center of different CNs. The atomic transition frequency 
is indicated by the dashed line in Fig. 13 (a) . 



center and near the wall inside, and at different distances 
outside achiral CNs of different radii. The relative density 
of photonic states (^{x) in Eq. ((13(1 was computed accord- 
ing to Eq. I(17() . The CN surface axial conductivity a^z 
appearing in Eq. I(17() was calculated in the relaxation- 
time approximation with the relaxation time 3xl0~^^ s; 
the spatial dispersion of ct^^ was neglected The free- 
space spontaneous decay rate was approximated by the 
expression ro{x)^a^2jox/h |35| . 

Figure[21(a) presents £,{x) for the atom in the center of 
the (5,5), (10,10) and (23,0) CNs. It is seen to decrease 
with increasing the CN radius, representing the decrease 
of the atom-field coupling strength as the atom moves 
away from the CN surface |23]- To calculate ]Cu(t)]^ in 
this particular case, we have fixed xa = 0.45 (indicated by 
the vertical dashed line), firstly, because this transition 
is located within the visible light range 0.305 < a; < 0.574, 
secondly, because this is the approximate peak position 
of ^(a;) for all the three CNs. The functions ](7„(t)]^ cal- 
culated are shown in Figure |21 (b) in comparison with 
those obtained in the Markovian approximation yield- 
ing the exponential decay. The actual spontaneous decay 
dynamics is clearly seen to be non-exponential. For the 
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FIG. 3: (Color online) (a) Fragment of the relative density of 
photonic states for the atom inside the (10,10) CN at distance 
of 3 A from the wall (the situation observed experimentally 
for Cs in Ref. |T3|'l. (b) The upper-level spontaneous decay 
dynamics for different atomic transition frequencies [indicated 
by the dashed lines in Fig.|3(a)] in this particular case. 



small radius (5,5) CN, Rabi oscillations are observed, in- 
dicating a strong atom-field coupling regime related to 
strong non-Markovian memory effects. Eq. (|21l) is satis- 
fied in this case. With increasing the CN radius, as the 
value of ^(0.45) decreases, Eq. (|^ ceases to be valid and 
the decay dynamics approaches the exponential one with 
the decay rate enhanced by several orders of magnitude 
compared with that in free space. Note that, though the 
distance from the atom to the CN surface is larger for 
the (23,0) CN than for the (10,10) CN, the deviation of 
the actual decay dynamics from the exponential one is 
larger for the (23,0) CN. This is an obvious consequence 
of the influence of a small neighboring peak in the (23,0) 
CN photonic density of states [Figure [3 (a)]. 

In Ref. formation of Cs-encapsulating single- wall 
CNs was reported. In a particular case of the (10,10) 
CN, the stable Cs atom/ion position was observed to 
be at distance of 3 A from the wall. We have simulated 
the spontaneous decay dynamics for a number of typical 
atomic transition frequencies for this case. Figure |3 (a) 
shows the density of photonic states and the five spe- 



cific transition frequencies xa (dashed lines) for which 
the functions |C„(i)|^ presented in FigureEKb) were cal- 
culated. Rabi oscillations are clearly seen to occur in the 
vicinity of the highest peak {xa ~ 0.22) of the photonic 
density of states. Important is that they persist for large 
enough detuning values xa ~ 0.21-^0.25. For xa — 0.30, 
the density of photonic states has a dip, and the decay 
dynamics exhibits no Rabi oscillations, being strongly 
non-exponential nevertheless. For xa = 0.45, the inten- 
sity of the peak of the photonic density of states is not 
large enough and the peak is too broad, so that strong 
atom-field coupling condition (|21|l is not satisfied and the 
decay dynamics is close to the exponential one. 

Figure ^ (a) shows the density of photonic states for 
the atom outside the (9,0) CN at the different distances 
from its surface. The vertical dashed lines indicate the 
atomic transitions for which the functions |Cu(t)|^ in 
Figures 21 (b), (c), and (d) were calculated. The tran- 
sitions Xa = 0.33 and 0.58 are the positions of sharp 
peaks (at least for the shortest atom-surface distance), 
while Xa = 0.52 is the position of a dip of the func- 
tion £,{x). Very clear underdamped Rabi oscillations are 
seen for the shortest atom-surface distance at a:A = 0.33 
[Figure0](b)], indicating strong atom-field coupfing with 
strong non-Markovity. For xa = 0.58 [Figure ^d)], the 
value of ^(0.58) is not large enough for strong atom-field 
coupling to occur, so that Eq. (|21|l is not fulfilled. As a 
consequence, the decay dynamics, being strongly non- 
exponential in general, exhibits no clear Rabi oscilla- 
tions. For 2:^1 = 0.52 [Figure^(c)], though ^(0.52) is com- 
paratively small, the spontaneous decay dynamics is still 
non-exponential, approaching the exponential one only 
when the atom is far enough from the CN surface. 

The reason for non-exponential spontaneous decay dy- 
namics in all the cases considered is similar to that taking 
place in photonic crystals 0. When the atom is closed 
enough to the CN surface, an absolute value of the rela- 
tive density of photonic states is large and its frequency 
variation in the neighborhood of a specific atomic tran- 
sition frequency essentially influences the time behavior 
of the kernel (|13|l of evolution equation (|12|l . Physically, 
this means that the correlation time of the electromag- 
netic vacuum is not negligible on the time scale of the 
evolution of the atomic system, so that atomic motion 
memory effects are important and the Markovian approx- 
imation in the kernel H13|) is inapplicable. 



V. CONCLUSIONS 

The effects we predict will yield an additional structure 
in optical absorbance/reflectance spectra (see, e.g., |2^ 
27]) of atomically doped CNs in the vicinity of the en- 
ergy of an atomic transition. Weak non-Markovity of the 
spontaneous decay (non-exponential decay with no Rabi 
oscillations) will cause an asymmetry of an optical spec- 
tral line-shape (similar to exciton optical absorbtion line- 
shape in quantum dots 36J). Strong non-Markovity of 
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FIG. 4: (Color online) (a) Relative density of photonic states for the atom located at different distances outside the (9,0) CN. 
(b, c, d) Upper-level spontaneous decay dynamics for the three atomic transition frequencies [indicated by the dashed lines in 
Fig. 2] (a)] at different atom-nanotube-surface distances. 



the spontaneous decay (non-exponential decay with fast 
Rabi oscillations) originates from strong atom-vacuum- 
field coupling with the upper state of the system split- 
ted into two "dressed" states. This will yield a two- 
component structure of optical spectra similar to that 
observed for excitonic and intersubband electronic tran- 
sitions in semiconductor quantum microcavities plll22l|. 

Summarizing, we have developed the quantum theory 
of the spontaneous decay of an excited two-level atom 
near a carbon nanotube. In describing atom-field inter- 
action, we followed an electromagnetic field quantization 
scheme developed for dispersing and absorbing media in 
Refs. 1^ 01 . This quantization formalism was adopted 
by us for a particular case of an atom near an infinitely 
long single- wall CN. We derived the evolution equation of 
the upper state of the system and, by solving it numer- 
ically, demonstrated a strictly non-exponential sponta- 
neous decay dynamics in the case where the atom is close 
enough to the CN surface. In certain cases, namely when 
the atom is close enough to the nanotube surface and 
the atomic transition frequency is in the vicinity of the 
resonance of the photonic density of states, the system 
exhibits vacuum-field Rabi oscillations, a principal sig- 



nature of strong atom-vacuum-field coupling, — a result 
already detected for quasi-2D excitonic and intersubband 
electronic transitions in semiconductor quantum micro- 
cavities (21j, 22^ and never reported for atomically doped 
CNs so far. This is the result of strong non-Markovian 
memory effects arising from the rapid frequency variation 
of the photonic density of states near the nanotube. The 
non-exponential decay dynamics gives place to the ex- 
ponential one if the atom moves away from the CN sur- 
face. Thus, the atom-vacuum-field coupling strength and 
the character of the spontaneous decay dynamics, respec- 
tively, may be controlled by changing the distance be- 
tween the atom and CN surface by means of a proper 
preparation of atomically doped CNs. This opens routes 
for new challenging nanophotonics applications of atom- 
ically doped CN systems as various sources of coherent 
light emitted by dopant atoms. 

Finally, we would like to emphasize a general char- 
acter of the results we obtained. We have shown that 
similar to semiconductor microcavities and photonic 
band-gap materials |^, carbon nanotubes may quali- 
tatively change the character of atom-electromagnetic- 
field interaction, yielding strong atom-field coupling — 
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an important phenomenon necessary, e.g., for quantum 
information processing js^, IH, |33- The present pa- 
per dealt with the simplest manifestation of this gen- 
eral phenomenon — vacuum-field Rabi oscillations in 
the atomic spontaneous decay dynamics near a single- 
wall carbon nanotube. However, similar manifestations 
of strong atom-field coupling may occur in many other 
atom-electromagnetic-field interaction processes in the 
presence of CNs, such as, e.g., dipole-dipole interac- 
tion between atoms by means of a vacuum photon ex- 
change cascade spontaneous transitions in three- 
level atomic systems |4ll |. 
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APPENDIX A: GREEN TENSOR OF 
A SINGLE-WALL CARBON NANOTUBE 

The form of the classical electromagnetic field Green 
tensor in Eqs. lO and H1U|) depends on the presence 
of external radiative sources (an atom in our case) and 
their position (inside or outside) with respect to a nan- 
otube. Since the problem has a cylindric symmetry, we 
assign the orthonormal cylindric basis {er,e^,ez} (see 
Figure ^ in such a way that is directed along the 
nanotube axis and ta = f = {rA, 0, 0}. 

We use the representation 

Gaz{r,rA,uj) = ^^VqVz + (5az^ g(r,rA, w), (Al) 

where a = {r^ip,z} and the function g(r, r^i, cj) is the 
Green function of the Hclmholtz equation. Substituting 
Eq. (HU into Eq. one straightforwardly obtains 



(A + fc2) g{r, rA,Lu) = -6{r - ta) 



with a known solution 



^ gifclr-r^l 

5o(r, rA,Lu) = — 

47r r-r^ 



(A2) 



(A3) 



satisfying the radiation condition at infinity (see, e.g., 
[s^). In our case, however, the functions G'az(r, r^, cj) 
and g(r, r^i, (jj) are imposed one more set of boundary 
conditions. They are the boundary conditions on the sur- 
face of the CN. Using simple relations (valid for r ^ ta 
under the Coulomb-gauge condition) 



(A4) 



and 



Hair,uj) = ^ eai3'r'^i3E^ir,u;) (A5) 

/3,7— r,(^,2: 



(eQ/37 is the totally antisymmetric unit tensor of rank 3), 
they can be derived from the classical electromagnetic 
field boundary conditions of the form 



E^ |r=_R„,+0 — E^ \r=R^„-0 = , 



(A6) 
(A7) 



47r 



H^\r=R^„+0 — H^\r=R^,^-0 — — CT zz{Rcn, ^^)E ^\r=R^ 



Ez\r = Rar.+0 - £lz\r = Rar.-^ 



H\r=f 







(A8) 
(A9) 



(spatial dispersion neglected) obtained in Ref. |23l |. 

Let rA>Rcn (the atom outside the CN) to be specific. 
Then, g{r, ta, w) may be represented in the form 

g{r r A to) - I 5o(r, r^, t^) + s^+Hr, ^) , r > R^n 
\ g'^^Hr^uj) , r < Ren 

(AlO) 

where go{'r,rA,'-^) is the point radiative atomic source 
function defined in Eq. (|A3|) and g^^\r,u;) are un- 
known nonsingular functions satisfying the homogeneous 
Helmholtz equation and the radiation conditions at in- 
finity. We seek them using integral decompositions over 
the modified cylindric Bessel functions Ip and Kp as fol- 
lows m 



and 



1 °° 

go{r,rA,u;) ^ ^ 



X / Ip{vr)Kp{vrA)e"^''dh, rA>r, 
Jc 



(A12) 



where Ap(h) and Bp(h) are unknown functions to be 
found from the boundary conditions (|A6|) - (jA9() in view 



of Eqs. IjAip . (|X4|) and HA5|I . v = v{h,uj) = VW^Ht^. 
The integration contour C goes along the real axis of the 
complex plane and envelopes the branch points ±fc from 
below and from above, respectively. 

The boundary conditions (|X6|l - (|X9|) with Eqs. IjAip . 
(|A4|) and (|A5I) taken into account yield the following 
two independent equations for the scalar Green func- 
tion (|ITO|) 

goir,rA,uj)\r=R,^ + g''^\r,u;)\r=R^^ = g^^\r,uj)\r=R,„, 



dg^~^^{r, L 
dr 



\r=Ra 



dg'- Hr,u}) 



dr 



(3{cu)(— + k']g(-\r,u;)\r=R^^^- 



dgo{r,rA ,uj) , 
dr 
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where /3(w) = 47ri (Tzz(i?cn, Substituting the inte- 

gral decompositions (|A11|I and (| A 12(1 into these equa- 
tions, one obtains the set of two simultaneous algebraic 
equations for the functions Ap{h) and Bp{h). The func- 
tion Ap{h) we need (we only need the Green function in 
the region where the atom is located) is found by solv- 
ing this set with the use of basic properties of cylindric 
Bessel functions (see, e.g., 43]). In so doing, one has 



ratio [see Eqs. lfTHl " lfT7|) ] 



where 



ImG'zz(rA,rA,^) 
ImG'^^(rA,rA,t^) ' 



(Bl) 



Ap{h) 



Rcnf3{Lo)v^l;{vRcn)Kp{vrA) 
(27r)2[l + /3(w) v^R„JpivR,^)Kp{vR,^)] 



These Ap{h), being substituted into Eq. (jAll|l . yield 
the function (7^+'(r, w) sought. The latter one, in view 
of Eq. IjAlOp . results in the scalar electromagnetic field 
Green function of the form 



g{r,rA,uj) = go{r,rA,uj) 



Rc 



(2^) 



I3{uj) v^ll{vR,,,)Kp{vrA)Kp{vr) 



Jc^ + V^RcnIp{vRcn)Kp{vR,n) 



Ahz 



dh , 



(A13) 

where > r > Ren- One may show in a similar way 
that the function g(r, r^,a;) for r<rA< Ren is obtained 
from Eq. IjAlSp by means of a simple symbol replacement 
Ip <-> Kp in the numerator of the integrand. 

Knowing g{r,rA,Lo), one may easily compute the 
components of the electromagnetic field Green tensor 
GQ2(r, r^,Li;) according to Eq. (|A1|I . 



Gzz(rA,rA,uj) 



1 g{r,rA,uj) 



(B2) 



according to Eq. (|A1|) . g{r,rA,i^) is given by Eq. (|A13|) 

for TA > r > Ren and by the same equation with the 
symbol replacement Ip Kp in the numerator of the 
integrand for r<rA< Ren ■ 



Substituting Eq. (|A13|I into Eq. HB2p . making differ- 
entiation and passing to the limit r ta, followed by 
the substitution of the result into Eq. (jBl|l . one has for 

rA > Ren 



3-Rr 



27rfc3 



X Im 



E 



(3iu)v^I^{vRen)K^ivrA) 

C H- Ren P{ll>) v'^ Ip{vRen)Kp{vRen) 



(B3) 



dh. 



APPENDIX B: DENSITY OF PHOTONIC 
STATES NEAR A SINGLE- WALL CARBON 
NANOTUBE 

The relative density of photonic states in an atomic 
spontaneous decay process near a CN is defined by the 



For TA < Ren , Eq. (|B3|I is modified either by the symbol 
replacement Ip Kp or by the Bessel function argument 
replacement <^Rcn in the numerator of the integrand. 
In dimensionless variables (|14|l . Eq. (|B3p is rewritten to 
give Eq. lfT7|l . 
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